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We propose and analyze a nanomechanical architecture where light is used to perform linear 
quantum operations on a set of many vibrational modes. Suitable amplitude modulation of a single 
laser beam is shown to generate squeezing, entanglement, and state-transfer between modes that are 
selected according to their mechanical oscillation frequency. Current optomechanical devices based 
on photonic crystals may provide a platform for realizing this scheme. 



The field of cavity optomechanics studies the inter- 
action between light and nanomechanical motion, with 
promising prospects in fundamental tests of quantum 
physics, ultrasensitive detection, and applications in 
quantum information processing (see [Ij for a review). 
One particularly promising platform consists of "optome- 
chanical crystals", with strongly localized optical and vi- 
brational modes implemented in a photonic crystal struc- 
ture [2J. So far, several interesting possibilities have been 
pointed out that would make use of multi-mode setups 
that can be designed on this basis. For example, suit- 
ably engineered setups may coherently convert phonons 
to photons [3j and collective nonlinear dynamics might 
be observed in optomechanical arrays [4J. Moreover, op- 
tomechanical systems in general have been demonstrated 
to allow writing quantum information from the light field 
into the long-lived mechanical modes [5-7J. The recent 
success in ground state laser-cooling [8j has now opened 
the door to coherent quantum dynamics in optomechan- 
ical systems. 

In this paper, we propose a general scheme for 
continuous-variable quantum state processing [9j utiliz- 
ing the vibrational modes of such structures. We will 
show how entanglement and state transfer operations can 
be applied selectively to pairs of modes, by suitable in- 
tensity modulation of a single incoming laser beam. We 
will discuss the limitations for entanglement generation 
and transfer fidelity, and show how to engineer the me- 
chanical frequency spectrum and pick suitable designs to 
address these challenges. 

Model. - We will first restrict our attention to a single 
optical mode coupled to many mechanical modes, such 
that the following standard optomechanical Hamiltonian 
describes the photon field a, the phonons hi of different 
localized vibrational modes (/ = 1, 2, . . . , A^), and their 
mutual coupling: 

H = -hAd'^d + ^^iW^i - ^^^^ So^ {bj^bi)^... . 

I I 

(1) 

Here we are working in a frame rotating at the laser 
frequency, with the detuning given by A = cjl — ^cav 
We omitted to explicitly write down the laser driving. 
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Figure 1: Schematic setup figure illustrating an optome- 
chanical crystal with localized vibrational modes, coupled by 
a common optical mode. The optical mode is driven by 
an amplitude-modulated laser beam to engineer frequency- 
selective entanglement, squeezing and state-transfer opera- 
tions of vibrations, as explained in the main text. 



and the coupling to the photon and phonon baths, with 
damping rates and F, respectively, although these will 
of course be taken care of in our treatment. The bare 
(single-photon) coupling constants ^ depend on the 
overlap between the optical and mechanical mode func- 
tions. They are generally on the order of co'cav^zPF/^, 
where L is an effective optical cavity length that reaches 
down to wavelength dimensions in photonic crystal cav- 
ities, and where xzpf = {h/2miVtiY^'^ is the mechanical 
zero-point amplitude of the respective mode (see Fig. [l] 
for the illustration of a setup). After going through the 
standard procedure of splitting off the coherent optical 
amplitude induced by the laser, d = a-\-Sd^ and omitting 
terms quadratic in Sd (valid for strong drive), we recover 
the linearized optomechanical coupling, 



H,nt = -h{Sd^Sd^)Y,9iibl^bi) =-h{Sd^Sd^)Y. (2) 



Here the dressed couplings gi 



.(0 



a can be tuned via 



the laser intensity, i.e. the circulating photon number: 
|<^| = V^phot 7 where we have taken a to be real- valued 
without loss of generality. We can now eliminate the 
driven cavity field (noting that 6d is in the ground state) 
by second-order perturbation theory. Provided we work 
at large detuning, | A| ^ 1]^, we retain a fully coherent, 
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light- induced mechanical coupling between all the modes, 



Lk 



(3) 



where Xi = bi-\-bl is the mechanical displacement in units 
of xzpF- Eq. ([3| may be viewed as a "collective optical 
spring" effect, coupling all the mechanical displacements. 
The couplings Jik = giQk/'^^ can be changed in-situ ei- 
ther via the laser intensity or the detuning. Note that 
if multiple optical modes are driven, the corresponding 
coupling constants will add. 

In general, the couplings in Eq. ([3| will induce quan- 
tum state transfer between mutually resonant mechanical 
modes, and entanglement at low temperatures (usually 
with the help of optomechanical laser cooling). These 
phenomena have been analyzed in a variety of schemes 
[TQHM] so far, typically with two mechanical modes of 
interest. 

General scheme. - However, here we have in mind 
a multi-mode situation for continuous variable quan- 
tum information processing. To this end, we are in- 
terested in having an efficient approach to selectively 
couple arbitrary pairs of modes, both for entanglement 
and state transfer. There are several desiderata to ad- 
dress for a suitable optomechanical architecture of that 
style : (i) The couplings should be switchable in a time- 
dependent manner; (ii) one should be able to easily se- 
lect pairs for operations; (iii) preferably, only one laser 
(or a limited number) should be involved; (iv) operation 
speeds should be large enough to overcome the effects of 
decay and decoherence; (v) one should be able to scale 
to a reasonably large number of modes. 

Static couplings as in Eq. ([3| could be used for selec- 
tive pairwise operations if one were able to shift locally 
the mechanical mode frequency, to bring into resonance 
only the two respective modes. In principle, this is doable 
via the optical spring effect, but would require local ad- 
dressing with independent laser beams. This could prove 
challenging in a micron-scale photonic crystal architec- 
ture, severely hampering scalability. 

Instead, we propose to employ frequency-selective op- 
erations, by modulating the laser intensity (and thus J) 
in a time-dependent fashion. Entanglement generation 
by parametric driving has been analyzed recently in var- 
ious contexts, including entanglement using supercon- 
ducting circuits [15j, trapped ions [16j, general studies 
of entanglement in sets of harmonic oscillators [T7H19J. 
optomechanical state transfer and entanglement between 
the motion of a trapped atom and a mechanical oscillator 
[20j and entanglement between mechanical and radiation 
modes [21j. Parametric driving can also lead to mechan- 
ical squeezing in optomechanical systems [22j. 

Let us consider two modes (1 and 2) for the mo- 
ment, where the coupling is 2hJ{t){Xi + ^2)^. Assum- 
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Figure 2: Squeezing and entanglement operations between 
three vibrational modes (Qi, Q2 and Q3) in an optomechanical 
setup, induced by harmonically modulating the laser intensity 
at a frequency uj. This schematic diagram shows the reso- 
nance regions for different processes starting to overlap when 
the effective light-induced coupling J is increased by chang- 
ing the laser modulation amplitude (Joveriap = (^^2 — ^i)/2 at 
first intersection of dashed lines). At large J the static part 
of the coupling also leads to unwanted swapping processes. 



ing a 100% amplitude modulated laser drive beam with 
|Qf(t)p = |amaxP cos^(a;t), wc have J{t) = Jcos^(cjt) = 
J[l + cos(2co't)]/2. The resulting time-dependent light- 
induced mechanical coupling can be broken down into 
several contributions, whose relative importance will be 
determined by the drive frequency u. The static terms, 
hJ{Xi + -^2)^7 will shift the oscillator frequencies by 
5Qj = 2J. In addition, they give rise to an off-resonant 
coupling (ineffective for — 1^21 ^ but with growing 
influence for \Qi — ^ J )• On the other hand, the 
oscillating terms contain 



(4) 



There are three important cases. A mechanical beam- 
splitter (state-transfer) interaction is selected for a laser 
drive modulation frequency uj = {Qi — 0.2) /2. After 
transforming the full Hamiltonian into the interaction 
picture with respect to Vti and 1^2, the resonant part 
then reads ^b.s. = hJib^^i + ^1^2)- In contrast, for 
UJ = {Qi + ^^2)72, we obtain a two-mode squeezing (non- 
degenerate parametric amplifier) Hamiltonian, i^ent = 
hj(bib2 -\-b\b2), which can lead to efficient entanglement 
between the modes. Finally, uj = selects the squeez- 
ing interaction for a given mode, Hsq = h{J/2){b^^ ^bf). 
These laser-tunable, frequency-selective mechanical in- 
teractions are the basic ingredients for the architecture 
that we will develop and analyze here. 

Limiting factors. - We now start to address the im- 
portant constraining factors limiting the fidelity of these 
operations, both for the two-mode and ultimately the 
multi-mode case. Full simulations incorporating all these 
effects will be discussed further below. At higher drive 
powers (as needed for fast operations), the frequency- 
time uncertainty implies that the different processes dis- 
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cussed above need not be resonant exactly any more, 
with an allowable spread \6uj\ < J. For example, the 
parametic instabilities occur for \uj — {fti -\- rtj)/2\ < J. 
At higher driving strengths, once these intervals start to 
overlap for different processes, selectivity is lost and the 
process fidelity suffers. On the other hand, at low opera- 
tion speeds quantum dissipation and thermal fluctuations 
will limit the fidelity. This dilemma is the essential prob- 
lem faced by a multi-mode setup, and we will discuss pos- 
sible schemes to address it further below. The schematic 
situation for three modes is illustrated in Fig. [2] 

In order to analyze quantitatively the full effects of de- 
coherence and dissipation, we employ a Lindblad master 
equation to evolve the joint state of the mechanical modes 
under the influence of light-induced time-dependent driv- 
ing. The evolution of any expectation value can be de- 
rived from the master equation and is governed by: 



ih 



[i,i^]) + ^(n, + l)^,(7^[6t]i 



^n,^,(7^[6,]i). 



dt 



Here H already contains the effective interaction ([3|, 
with modulated time-dependent couplings, Tj are the 
mechanical damping rates of the vibrational modes, and 
Uj their equilibrium occupations at the bulk tempera- 
ture. The dissipative term in the first line describes 
damping (spontaneous and induced emission), and the 
second line refers to absorption of thermal fluctuations. 
The relaxation superoperators are defined by 7^ [6^] A = 
Ab—b^bA/2—Ab^b/2 (in contrast to the equation for p). 
For the quadratic Hamiltonian studied here, the equa- 
tions for correlators, such as (b]b^ \ remain closed, and 



these (together with averages y^j^ ^ire sufficient to de- 
scribe the time-evolution of the Gaussian quantum states 
that will be produced in the course of the dynamics. 

In order to analyze selective entanglement, we evaluate 
the logarithmic negativity 



E^r{p) = log2 



(5) 



as a measure of entanglement for any two given modes {A 
and 5), where pab is the state of these two modes, and 
the partial transpose Ta acts on A only. For Gaussian 
states. Ex can be calculated by obtaining the symplectic 
eigenvalues of the covariance matrix of the two modes' 
positions and momenta [23J. 

In Fig. |3] we show the results of numerical simulations 
for a situation with three vibrational modes, two of which 
are to be entangled in the presence of the third one. 
The entanglement first grows and then saturates at later 
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Figure 3: Selective entanglement between two vibrational 
modes, in an array of three, produced by an amplitude- 
modulated laser beam driving a common optical mode. All 
modes are coupled to a thermal bath and assumed to be 
cooled to their ground state via side-band cooling prior to 
the entanglement operation. (a,b) Color plots showing an 
entanglement measure (the logarithmic negativity) between 
modes 1 and 3, as a function of the driving strength J and 
the modulation frequency cj, at a fixed time. In (a) the vibra- 
tional modes are quite close [Qi = (1, 1.025, 1.075)^^], which 
leads to overlapping resonance regions for different processes 
(as in Fig. [2]). In (b), the frequency spacing is increased by a 
factor of four, significantly increasing the selectivity, (c) De- 
pendence on driving strength, i.e. cut along lines indicated 
in (a,b), see main text, (d) Dependence on temperature and 
mechanical quality factor, (e) Same as (d) but without pre- 
cooling, i.e. starting from thermal equilibrium. [Contour lines 
in (a,b): entanglement between 1&2 or 2&3 larger than 1.0, 
blue arrows at Qj; Parameters for (a,b,c): n — 100, Q — 10^; 
(d,e): J/^ = 0.003; (b,d,e): n^ = (1, 1.1, 1.3)1^] 



times, while the phonon number continues to grow expo- 
nentially. The plots show the entanglement evaluated at 
a fixed late time (t = 5.6/ J), as a function of parameters. 
One clearly sees the features predicted above, i.e. the un- 
wanted overlap between different entanglement processes 
at higher driving strengths (Fig. [3|i,b). Increasing the 
vibrational frequency spacing suppresses these unwanted 
effects (Fig.|3]3). The dependence on the driving strength 
itself is displayed in more detail in Fig. [3]:. There, the 
threshold J = Vfi for entanglement generation at finite 
temperature is evident, as is the loss of entanglement at 
large J. Finally, Fig. [3]i,e shows the dependence on tem- 
perature and mechanical quality factor. It indicates that 
this scheme should be feasible for realistic experimental 
parameters (see below). Note that the light-induced dis- 
sipation [24J effectively adds to the intrinsic decoherence 
rate Tn the rate F^p^ ^ gla^n//S? = 2J(/^/A). This 
is suppressed by a factor tv/A which can in principle be 
made arbitrarily small for larger detuning (at the expense 



of higher circulating photon number o? to keep the same 
J). For the reahstic experimental parameters quoted be- 
low, we have /^/A = 1/80, such that we have been able 
to neglect the effects of F^p^. 

Larger arrays. - We now turn to the situation with 
an array of many modes. It is clear that having evenly 
spaced mechanical frequencies is impossible without tak- 
ing any further precautions. This is because then the 
state transfers between adjacent modes would all be ad- 
dressed at the same modulation frequency. In fact, there 
seems to be no layout that allows for selection of arbi- 
trary pairs, avoids resonance overlap, and does not re- 
quire a frequency interval that grows exponentially with 
the number N of modes. While one may still realize 
small arrays in this way, in the limit of large TV another 
approach is needed. 

The scheme (Fig. |4| that solves this challenge in- 
volves an auxiliary mode at ^laux, removed in frequency 
from the array of "memory" modes which now may have 
evenly spaced frequencies in an interval [^^min, ^max]- All 
the pairwise operations will take place between any se- 
lected memory mode and the auxiliary mode. Then, 
the state transfer resonances are in the band [(l^aux — 
^max)/2, (^aux— ^min)/2], and entanglement is addressed 
within [{Qmin + ^aux)/2, (l^max + l^aux)/2]. To make this 
work, one needs to fulfill the mild constraint 21]max — 
^min < ^aux < Sl^min- State transfer between two mem- 
ory modes now is performed in three steps (swapping 
1-aux, aux-2 and aux — 1), as is entanglement (swap 
1-aux, entangle aux-2 and swap aux-1). Note that this 
overhead does not grow with the number of memory 
modes. Fig. [5] shows the state transfer between an auxil- 
iary mode (originally prepared in a squeezed state) and 
one of the memory modes. This scheme can potentially 
be expanded, with the auxiliary modes grouped into ar- 
rays (Fig. [4]:), and several of such 2D blocks could again 
be connected via further "higher-order" auxiliary modes, 
in a hierarchical fasion. 

Implementation. - Regarding the experimental imple- 
mentation, in principle any optomechanical system with 
several long-lived mechanical modes can be used as a 
starting point. One promising platform is based on pho- 
tonic crystals ("optomechanical crystals"), as introduced 
by Painter et al. [2j. These would be very well suited 
for the scheme presented here, due to their design flexi- 
bility, particularly of two-dimensional structures, and the 
all-integrated approach, as well as the very large optome- 
chanical coupling strength. Given the current coupling 
strength achieved there ^o/27r ~ IMHz, as well as 
a detuning of A/Q = 10 and around 2000 photons cir- 
culating inside the cavity (a number reached in recent 
experiments), we can estimate the induced coupling to 
approach the damping rate, J ~ F. This corresponds 
to the threshold for coherent operations, provided one 
were to cool down the bath to /c^^bath < This is in 
principle doable (at 20mK), but will likely run into prac- 
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Figure 4: (a) Layout for an optomechanical array with an 
auxiliary vibrational mode, circumventing the problem of fre- 
quency crowding (see main text), (b) Corresponding mechan- 
ical frequency spectrum. The laser modulation frequency has 
to lie within certain intervals to select entanglement or state- 
transfer operations between the "memory modes" and the aux- 
iliary mode, (c) Extension to the second dimension with op- 
erations "around the corner". 
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Figure 5: State transfer of a squeezed state (Ax^ 
from the auxiliary to the second memory mode (all in the 
ground state due to prior side-band cooling) (a) Dependence 
of the transfer fidelity on driving strength and pulse time, 
with a swap pulse ideally occuring at Jt = 7r/2. For long 
pulse times dissipation hampers the transfer, for large driving 
strengths resonances overlap, medium times are optimal, (b) 
Wigner density of the initial squeezed state in the auxiliary 
mode, and (c) after the transfer, (d) Maximization of transfer 
fidelity over parameter ranges used in (a). Parameters: Qj = 
(1.0, 1.1, 1.2, 2.0)Q, n = 100, Q = 10^ J/Q = 0.007 (b,c). 



tical difficulties due to the re-heating of the structure 
via spurious photon absorption or other effects. Oth- 
erwise, at finite bath temperatures corresponding to a 
thermal occupation n ~ /^B^bath/^^, the light intensity 
must be increased by a factor n, towards J ^ Fn, to 
speed up operations and thereby fight thermal decoher- 
ence. In that case, the vibrational ground state would 
be prepared at the start of the pulse sequence via laser- 
cooling, as demonstrated in [8j. 

In these devices, several localized vibrational and op- 
tical modes can be produced at engineered defects in an 
otherwise periodic array of holes cut into a free-standing 
substrate (e.g., made of silicon). Evanescent optical and 
vibrational waves connect adjacent modes (via photon 



5 



and phonon tunneling, respectively). The typical photon 
tunnel coupling for modes spaced apart by several lat- 
tice constants is [4J in the range of several THz. Thus, 
hybridized optical modes will form, one of which can be 
selected via the laser driving frequency as the active com- 
mon optical mode (the others remaining idle). At the 
same time, the vibrational modes' frequencies can either 
be designed to be different or to be equal, in which case 
delocalized hybridized mechanical modes are produced. 
Recently it was shown that a 'snowflake' crystal made 
of connected triangles (honeycomb lattice) possesses a 
simultaneous photonic and phononic (pseudo-)bandgap 
and thus supports wave guides (line defects) and local- 
ized defect modes with optomechanical interaction [3j. 
Placing point defects (heavier triangles/thicker bridges) 
in the middle of such a crystal structure, a tight binding 
analysis indicates that the desired mechanical frequency 
spectrum (Fig. [4| can be generated in principle. 

In any given system there will be limits to the de- 
sign of the mechanical spectrum. We briefly mention an- 
other option for improving the operational fidelity: pulse 
shaping and optimal control. Essentially, one wants to 
make sure that the Fourier transform of the coupling J{t) 
(or, equivalently, of the time-dependent laser intensity), 
does not contain spectral weight at any of the resonances 
00^ = {fti ± %)/2, except for the selected one. This im- 
plies that the pulse duration is larger than the inverse of 
the smallest spacing of such resonances. Optimal control 
techniques (as in [17j) could be employed to numerically 
search for the optimal pulse shape. 

Finally, one essential ingredient of any such architec- 
ture will be read-out. Some time ago, we have pointed 
out [25] how to produce a quantum-non-demolition read- 
out of the quadratures of mechanical motion in an op- 
tomechanical setup. A laser beam impinging onto the op- 
tical resonance (detuning A = 0) is amplitude-modulated 
at the mechanical frequency Qj of one of the modes. The 
reflected light carries information only about one quadra- 
ture e'^'-^bj-\-e~'^'-^by Its phase (p is selected by the phase of 
the amplitude-modulation, while the measurement back- 
action perturbs solely the other quadrature. In that way, 
by repeated measurements all the joint correlators of po- 
sitions and momenta of the mechanical modes may be 
read out, e.g. in order to verify the fidelity of the op- 
erations discussed above. If desired, taking measure- 
ment statistics for continuously varied quadrature phases 
would also allow to do full quantum-state tomography of 
the set of vibrational modes, and thereby ultimately pro- 
cess tomography. 

Conclusions. - The scheme described here would en- 
able coherent scalable nanomechanical state processing in 
optomechanical arrays. It can form the basis for generat- 
ing arbitrary entangled mechanical Gaussian multi-mode 
states. An interesting application would be to investigate 
the decoherence of such states due to the correlated quan- 
tum noise acting on the nanomechanical modes. More- 



over, recent experiments have shown in principle how ar- 
bitrary states can be written from the light field into the 
mechanics jSHZ]- These could then be manipulated by 
the interactions described here. Alternatively, for very 
strong coupling go > hi, non- Gaussian mechanical states 
[26] could be produced, and the induced nonlinear in- 
teractions (see e.g. [27l [28]) could potentially open the 
door to universal quantum computation with continuous 
variables [H] in these systems. 
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